Ripple edge engineering of graphene nanoribbons 
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It is now possible to produce graphene nanoribbons (GNRs) with atomically defined widths. GNRs 
offer many opportunities for electronic devices and composites, if it is possible to establish the link 
between edge structure and functionalisation, and resultant GNR properties. Switching hydrogen 
edge termination to larger more complex functional groups such as hydroxyls or thiols induces strain 
at the ribbon edge. However we show that this strain is then relieved via the formation of static out- 
of-plane ripples. The resultant ribbons have a significantly reduced Young's Modulus which varies 
as a function of ribbon width, modified band gaps, as well as heterogeneous chemical reactivity 
along the edge. Rather than being the exception, such static edge ripples are likely on the majority 
of functionalized graphene ribbon edges. 

PACS numbers: 81.05.ue, 62.20.de, 81.40.Jj, 31.15.E-, 81.05.uj 
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I. INTRODUCTION 

Since the start of the graphene boom in 2004 1 , the 
simplistic picture of a perfect flat carbon monolayer has 
been refined. Notably dynamic rippling of graphene has 
been demonstrated via electron microscopy 2 - and molec- 
ular dynamics calculations^, and there have been first 
discussions of the possibility for stable static ripples in 
graphene^ 5 -. In reality graphene is not an infinite plane 
but is constrained by edges. Graphene nanoribbons 
(GNRs) have different properties from the infinite bulk 
material, notably they can display a finite band gap as 
a function of ribbon width^. In addition the physical 
and chemical behaviour of the one-dimensional edges will 
be superimposed on that of the bulk two-dimensional 
graphene. Graphene nanoribbons can be produced by 
unzipping carbon nanotubes 7 , lithography 8 ? 9 , etchingiil 
and controlled chemical bottum-up methods 11 , and 
offer great potential both for nanoelectronics 12 and 
nanocompositesii Edges are an easily accessible way 
to chemically functionalise the graphene and hence mod- 
ify its properties. We show here that careful design of 
graphene edges allows us to define the graphene nanorib- 
bon properties. 

The simplest way to saturate GNR edge dangling 
bonds is via hydrogen termination, which is the standard 
approach in GNR modellin g 14 ? 15 . However other termi- 
nating heteroatoms can be imagined (e.g. N, O), and an 
intriguing example of this was a recent theoretical study 
which found F-terminated armchair GNRs to be more 
stable when twisted helically^. 

However to develop a more realistic picture of the 
possibilities of GNR edge chemistry, more complex 
termination groups have to be investigated. A good 
example are hydroxyl groups (OH), which as well as 
being bulkier than simple heteroatoms also show more 
complex chemical interaction between themselves. Such 
groups have been proposed as a way to introduce strain 



along the graphene edge 17 , in order to tune GNR 
electronic properties such as bandgap 18 ? 19 . However 
two dimensional layered materials such as graphene 
have alternative mechanisms for relieving edge-induced 
strain, namely structural deformation into the third 
dimension via rippling or buckling 5 . This additional 
degree of freedom adds significant richness to graphene 
edge chemistry, which we investigate here for infinitely 
long armchair graphene nanoribbons (AGNRs) using the 
example of hydroxyl functionalisation. In the current 
study we show how -OH termination of different width 
AGNRs modifies their structure, strain and stability. 
We then demonstrate the effect of such functionalisation 
on band gap, chemical reactivity to metal deposition, 
and Young's Modulus, and generalise to other functional 
groups and ribbon types. This study points the way 
towards "edge termination engineering" as a way to 
create GNRs with custom designed properties. 



II. METHOD 

We performed density functional theory calculations 
under the local density approximation as implemented 
in the AIMPRO cod o 20 ? 21 . The calculations were carried 
out using super cells, fitting the charge density to plane 
waves within an energy cut-off of 200 Ha. Electronic level 
occupation was obtained using a Fermi occupation func- 
tion with kT = 0.04 eV. Relativistic pseudo-potentials 
are generated using the Hartwingster-Goedecker-Hiitter 
scheme 22 . These functions are labelled by multiple or- 
bital symbols, where each symbol represents a Gaussian 
function multiplied by polynomial functions including all 
angular momenta up to maxima p (I = 0, 1) and d (I = 
0, 1, 2). Following this nomenclature, the basis sets used 
for each atom type were pdddp (C), ppp (H) and dddd 
(O), resulting in 38 independent functions for carbon, 12 
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FIG. 1. (Color online) H-terminated AGNR, unit cells shown 
by dotted frames, definition of armchair ribbon width shown 
by numbers 24 , width 6. C atoms are shown grey and H atoms 
are white. 
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for hydrogen and 40 for oxygen. A more detailed account 
of the basis functions can be found elsewhere 23 . A Bloch 
sum of these functions is performed over the lattice vec- 
tors to satisfy the periodic boundary conditions of the 
supercell. 

Supercell sizes have been checked and chosen to be 
sufficiently large (y- and z-distance between ribbons > 12 
A) to avoid interaction with neighbouring GNRs. A fine 
k-point grid was chosen of the form 12/(n-ao)xlxl with 
n G N where ao is the length of a fundamental unit cell 
along the ribbon axis in the supercell (see FigJT]), which 
gives energies converged to better than 10 -5 Ha. For 
rippled GNRs the unitcell was doubled along the ribbon 
(2 • ao) to satisfy the periodic conditions of the supercells 
(longer period ripples using n • ao cells, n > 2 were also 
tested but found to be less stable) . For all structures 
the atom positions and the lattice parameters have been 
fully relaxed. The definition of ribbon widths of AGNRs 
is given by Cervantes- Sodi et al^ and is used in the text 
and Fig J3] and Fig |6] (definition see FigJT]). 





III. RESULTS AND DISCUSSION 

Larger functional groups on the GNR edge can intro- 
duce a range of possible inter-group interactions, includ- 
ing steric hindrance, Coulombic repulsion, dipole-dipole 
interactions and hydrogen bonding. Taking the case of 
hydroxyl groups we attempt to address these systemati- 
cally. We have modelled fifteen different structural pos- 
sibilities but focus here on three structures of hydroxy- 
lated edges for GNRs with widths from 4 to 20 (refered 
to hereafter as structures A, B and C). The relevant re- 
laxed edge structures are shown in Fig|2] and associated 
structural parameters in Table ID 

We start with structure A, symmetry constrained to lie 
in the plane with symmetrically paired -OH groups, re- 
moving any possibility for hydrogen bonding. Calculated 
bond lengths are dc-o = 1.33 A and do-H = 0.96 A. 
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FIG. 2. (Color online) (a)-(c): three OH-terminated arm- 
chair edge structures: A (metastable base symmetric struc- 
ture), B (metastable with hydrogen bonding) and C (stable 
ground state structure with static rippling), (d)-(e): stable 
OH-terminated AGNR of type C (width 6), (d) perspective 
view of ribbon and (e) side view with clearly visible rippled 
edge. Dotted lines are a guide to the eye for the GNR. u 
stands for "up" and d for "down" for the rippled edge. C 
atoms are in grey, O atoms in red and H atoms are white 
(same for Figgl FigEJand FigE}. 
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FIG. 3. (Color online) (a): Energy difference AU between structure A (black) and structure B (green) referring to structure 
C (blue) versus ribbon width from 4 to 20. (b): Strain along the ribbon (x-axis see Fig[T]), strain is referenced to a perfect flat 
graphene sheet (e = %), H-termination (red), A (black), B (green) and C (blue). Fits with e(x) — -% + | + c (fit parameter 
in Supplementary Materials^). The top scale of (a) and (b) in A is taken from perfect flat graphene with C-C bondlengths of 
d c -c — 1.41 A as guide for the eye (also for Fig|6]). 



Structure 


dc-o (A) 


do-H (A) 


do-H (A) 


A (flat) 


1.33 


0.96 


- 


B (flat) width 4 
width 20 


1.36 
1.35 


1.05 
1.09 


1.26/1.41 
1.11/1.14 


C (rippled) width 4 
width 20 


1.36 
1.36 


1.02 
1.02 


1.40/1.63 
1.45/1.69 


Phenol (C 6 H 6 0) 


1.35 


0.97 


- 



TABLE I. Overview of calculated bond lengths of the AGNRs 
from width 4 — 20 for structures A, B and C. O- • H bond 
lengths for Structure B and C vary with ribbon width between 
the limits given in the table. Structure C is the stable ground 
state. Calculated Phenol values included for comparison. 



This structure was recently discussed as a possible con- 
figuration for inducing strain in AGNRsii. Indeed we 
find that it induces strains e of up to 9 % for small AG- 
NRs (see Figj3jb). However this structure is extremely 
unstable. Strain along the ribbon has been calculated as 
follows (see also FigJT]): 

^0 ^Graphene 



(^Graphene 



with aQraphene = 4.23 A our optimized DFT (LDA) lat- 
tice parameter along the armchair direction of free stand- 
ing graphene, in good agreement with experimental val- 
ues of w 4.2 A2£. 

We next introduce hydrogen bonding between the hy- 
droxyl groups by breaking the in-plane symmetry, giving 
the most stable planar configuration, structure B. Here 
dc-o = 1.36 A with hydrogen forming one strong co- 
valent bond do-H — 1-05 — 1.09 A and one weaker hy- 
drogen bond do-H — 1-H — 1-41 A. This new hydrogen 
bond lowers the system energy by typically 3.7 eV/unit 



cell compared to Structure A, and reduces slightly the in- 
duced strain (still up to 6.7% for small AGNR compared 
to perfect graphene, as can be seen in Figj3jb). 

However by additionally breaking the planar symme- 
try we reach the energetically most stable configuration 
for hydroxyl terminated AGNRs, structure C. The -OH 
groups displace out of plane pairwise, creating a static 
sinosoidal ripple along the AGNR edge (see Figj2jd,e). 
By displacing out of plane the structure releases up to 
1.46 eV/unitcell (Figj3ja) and the strain in the ribbon 
is relieved (Fig|3jb), returning to values similar to hy- 
drogen terminated AGNRs. The carbon pairs for each 
hexagonal ring are displaced up and down alternately 
(longer period oscillations tested in supercells with n • ao 
and n > 2 at a variety of ribbon widths were found to 
be less stable). While dc-o stays largely unchanged at 
1.36 A, the covalent do-H extends to 1.02 A, while the 
rippling allows the hydrogen bond lengths to increase to 
do-H — 1-40 — 1.69 A. These lengths are approaching 
our calculated bondlengths for phenol (dc-o — 1-35 A, 
do-H = 0.97 A, Table P. 

The out of plane deformation mode can be understood 
as an elastic response to a ID edge line tension applied to 
a rigid 2D graphene sheet (our calculated Young's Modu- 
lus for graphene is 1.08 TPa), where it is clearly energet- 
ically favourable to buckle the edge rather than stretch 
the whole sheet. However non-planar sheet deformation 
decreases the 7r-orbital overlap, and this enthalpic driv- 
ing force of graphene to remain fiat localises the rippling 
along the ribbon edge. This can be seen in FigH] where 
the ripple amplitude increases from 1.63/2 = 0.82 A 
(width = 4) to 1.78/2 = 0.89 A for large ribbons (width 
> 12), and where out-of-plane displacements remain lo- 
calised within w 2.3 A of the graphene ribbon edge. The 
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FIG. 4. (Color online) Front view of -OH terminated AGNRs 
(stable structure C), (a) width 6, (b) width 12 and (c) width 
18. 
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ripple amplitude is inversally proportional to the AGNR 
strain and only indirectly coupled to the width of the rib- 
bon, consistent with a picture of a constant buckled edge 
length more able to dilate the graphene when its basal 
plane area is smaller. 

No energy barrier was found between the rippled C and 
flat B structures, demonstrating that configuration B is a 
metastable maximum (see Supplementary Materials^). 
The energy barrier calculations can be summed to es- 
timate the ripple inversion barrier, giving a minimum 
barrier of AUsarrier = 0.59 eV/ao to invert one edge 
wave. This barrier is too high for thermal activation at 
room temperature, and structure C can be considered as 
a static edge rippling. 

Thus we find that hydroxylated AGNRs will show 
spontaneous static edge ripples, due to an edge-strain 
induced out-of-plane deformation mode. As a result 
the strain induced by the rippled hydroxylated edge 
is for wide ribbons less then 0.35 % (width > 20), 
and very close to that of H-terminated AGNRs of the 
same width (< 0.18 %). We note that we found no 
significant interaction between the two edges of different 
OH-terminated AGNRs (as reflected in the total energy 
or strain), suggesting they are largely decoupled. 

We next investigate how static ripples modify the rib- 
bon properties, comparing structures B and C to deter- 
mine the influence of rippling, and C and H-terminated 
ribbons to determine the influence of the hydroxyl func- 
tional group. Starting with electronic properties, we 
obtain good agreement between our calculated LDA 
band gap for flat H-terminated AGNRs and previous 
literature^. Bandgap is inversely proportional to ribbon 
width, superimposed with strong 3N periodicity. This 
alternation is explainable via Clar sextet theory^ and 
Fermi wavelength^. The calculated band gap for struc- 
tures B and C (see FigjSJa) are almost superposed, show- 
ing that edge rippling does not appear to affect the gap 
for -OH termination, which is instead dominated by the 
choice of edge functional group. Once again we observe 



FIG. 5. (Color online) Pd atom addition to a -OH terminated 
AGNR (width 10). (a) "Top-ridge" site, (b) "Valley" site, 
which is AU = —0.4 eV more stable. 



3N periodicity^ but the difference with H-termination is 
not a simple phase shift of the periodicity. For 3N + 2 
(N = 1,2,..) widths the bandgap is similar to the H- 
terminated case, but for 3N and 37V + 1 the difference 
is significant, with bandgap fluctuations of up to 50%. 
Thus ribbon band gap appears highly sensitive to choice 
of edge functional group, however the gap is not sensitive 
to out-of-plane edge rippling for -OH termination. 

Rippling may be expected to modify the chemical re- 
activity of the ribbon surface, and to estimate this we 
calculated possible bonding sites for single Pd atom ad- 
dition on a AGNR (width 10). Pd sits above C-C bond 
centre o 29 i 30 , confirmed in our calculations. In general 
binding nearer the ribbon edge is more stable. For rip- 
pled edges we find Pd atoms at the ribbon edge are 0.4 
eV more stable in a concave "valley" site of an edge ripple 
than on the convex "ridge top" (see FigJSJ). Indeed the 
concave "top" site is even 0.1 eV less stable than a flat 
site at the ribbon centre. We expect that with enough 
activation energy the Pd atoms can migrate to the edge 
and will sit in the valleys of the rippled edge structure. 
Thus the ripple-induced changes in surface curvature pe- 
riodically modify the surface reactivity of the graphene, 
suggesting interesting changes in absorption behaviour 
for chemisorbed and physisorbed species, and potentially 
important geometric effects for metal contact deposition. 

The largest change due to edge rippling we find is in 
AGNR mechanical behaviour. We determined Young's 
Modulus (E) of AGNRs along their length using the ap- 
proach described by Zeinalipour-Yazdi et al 31 . Induc- 
ing strain up to ±3 % along the ribbon the total ener- 
gies have been calculated (7-points). Energy difference 
(AU) vs. strain (e) was fitted with a quadratic function 
AU(e) = ae 2 . E is given by E = ^, where V = a Q xwxh 
is the volume of the ribbon section composed of length 
ao, width w taken as the H-H distance between the edges, 
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FIG. 6. (Color online) (a): Energy gap versus width for different terminated AGNR. H-terminated AGNR (red), structure B 
(green) and structure C (blue), (b): Young's Modulus E along the AGNRs (in case of graphene along the armchair direction). 
E fitted with ^=f + -r + cas guide for the eye (fit parameter in Supplementary Materials^). 



and h = 3.35 A. We find E for flat infinite graphene to be 
1.08 TPa, in good agreement with literature values (1.09 
TPa^i, 1.05 TPa22). However finite width H-terminated 
and edge-rippled OH-terminated AGNRs show signifi- 
cantly smaller Youngs' Modulus (see FigjSJb). In both 
cases modulus is approximately inversely proportional to 
ribbon width, with some non-linearity at smaller widths 
when edge effects start to dominate. Both edge termina- 
tions extrapolate to the ideal graphene value at infinite 
width (H-terminated: 1.05 TPa, OH-rippled: 1.03 TPa), 
and as expected the modulus of -H terminated AGNRs 
drops towards the value for cis-polyacetylene at small 
widths. However surprisingly, changing the edge func- 
tionalisation from -H to rippled -OH groups significantly 
decreases the Youngs' Modulus of the ribbon. This re- 
duction is a direct result of the rippling (indeed calcula- 
tions for the unstable flat -OH ribbons structure B, actu- 
ally show a slight increase in Youngs' Modulus over the 
-H terminated case, since in this case the ribbon is under 
slight tension (see Supplementary Materials^)). We note 
that as for the band gap, a periodicity of 37V can be seen 
in the values of E reflecting variations in C-C bonding 
with width. These changes in E as a function of both 
ribbon width, and edge termination, will be critical in 
graphene nanocomposite design. Modifying ribbon edge 
functionalisation is an obvious way to chemically bind 
nanoribbons into a host polymer matrix, but these re- 
sults show that changing the functionalisation groups can 
decrease the Young's Modulus of the nanoribbon itself as 
much as 40 %. 

In order to determine how general the static edge rip- 
pling effect is, we next examined other chiralities. Hy- 
droxylated zig-zag GNRs (ZGNRs) do not exhibit out- 
of-plane edge rippling, due to the larger spacing between 
edge C atoms (2.44 A) which results in sufficient spacing 
between the -OH groups to keep edge strain below ±0.2 
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FIG. 7. (Color online) Chiral -OH terminated GNR. (a) per- 
spective view, (b) side view, (c) front view. 
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FIG. 8. (Color online) AGNR (width 7) -SH terminated, (a) 
perspective view, (b) side view, (c) front view. In grey C 
atoms are pictured, H atoms in white and S atoms in yellow. 
u stands for "up" and d for "down" for the rippled edge. 



%. Indeed, in principle edge tension might be expected 
to induce rippling within the ribbon centre, the inverse 
effect of AGNRs, however the rigidity of the graphene 
sheet precludes this. Note however this does not preclude 
rippling in ZGNRs for larger functional groups. Calcula- 
tions for hydroxylated chiral edges show an intermediate 
effect, with localised edge rippling around armchair-like 
sections which rapidly decays in zig-zag sections of the 
edge (see FigJT]). Thus these results show that the major- 
ity of GNR edge types, when hydroxylated, will exhibit 
static rippling. 

We also examined the dependence of static edge rip- 
pling in AGNRs on functional group, replacing -OH with 
-F, -CI, and -SH. In all cases the ribbons exhibited the 
same periodic edge rippling as the hydroxylated edges 
(see Fig J8]and Supplementary Materials 25 ). Thus a key 
finding of this study is that flat GNR edges, as exem- 
plified by hydrogen termination, appears to be the ex- 



ception rather than the rule. The majority of graphene 
nanoribbon chiralities and functionalisations we have ex- 
amined undergo spontaneous out-of-plane static rippling. 



IV. CONCLUSIONS 

This fundamental study provides a first picture of un- 
derlying physics for adding complex functional groups 
to graphene edges and in particular AGNRs. Hydroxyl 
groups do not induce large strain in GNRs, since the 
strain is compensated by static ripple formation along 
the ribbon edge. Such ripples form in the majority of 
ribbon chiralities and functional groups we have exam- 
ined, suggesting flat ribbon edges as observed for hydro- 
genation may be the exception. This also seems to rule 
out edge functional groups as a simple way of inducing 
large strains in the graphene basal plane. For the first 
time we calculate Young's Modulus for infinite AGNRs 
with -H and -OH terminations, and show that edge rip- 
ples can drastically modify the mechanical properties and 
chemical reactivity of the ribbon. We believe this could 
have a big impact on engineering devices and composites 
with embedded GNRs. We find that the band gap is not 
sensitive to static edge rippling for -OH termination, but 
can change by up to 50 % depending on choice of edge 
functional group. 

There are various ways to match a ID line tension 
along a ribbon edge against a 2D surface strain in the 
basal plane. In the case of graphene, with a Young's 
Modulus of 1.08 TPa, the in-plane resistance to ten- 
sion is high and edge compression is the energetically 
favoured solution, resulting in static ripple formation at 
the edge. At the same time there is an energetic cost asso- 
ciated with rippling induced disruption of the graphene 
7r-network, and hence the ripples remain localised near 
to the ribbon edge, reducing the effective flat basal plane 
width of the ribbon. In the case of chiral and aperiodic 
edges, static edge rippling is likely to lead to weak locali- 
sation and a decrease in the intravalley scattering length, 
consistent with experimental observations in graphene 
flakes^. 

The discussion thus far concerns free-standing 
graphene. On substrates the modulating height of the 
hydroxylated edge ripples will result in periodic mod- 
ulation in the graphene-substrate spacing and hence 
the substrate-induced potential felt by the nanoribbon, 
increasing electronic localisation effects. Variation in 
graphene-substrate spacing at the ribbon edge could also 
facilitate impurity intercalation beneath the graphene. 
Finally we note that this static edge rippling behaviour 
is a fundamental response of a 2D layered system to ID 
edge strain, and as such is also likely to also be of im- 
portance in the range of new monolayer materials under 
development such as BN, M0S2 and NiTe^. Strain com- 
pensation through rippling is also likely to be important 
at other graphene interfaces such as grain boundaries. 
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